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^vq We have measured second- and third-order cumulants in UAl data (pp collisions 

at 630 GeV/c). Rather than quoting numerical values for source parameters, we 
have used these in three checks to test the "quantum statistics" theory for consis- 
■ tency over these cumulants. In the process, we have found a method for folding 

theoretical correlation functions with experimental one-particle distributions. Our 
O preliminary results appear to indicate that, for the specific tests performed, the data 

ft, ■ contradicts the theory. 
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1 Introduction 

Pion interferometry has been a part of particle physics for several decades [0]. The main 
sub-branch of this science, concerning itself with Bose-Einstein correlations, endeavours 
to elicit information on the size, shape and temporal evolution of the source emitting 
pions. This is based on an analogy between optical intensity interferometry and quantum 
mechanical interference between incoherent pion amplitudes. 

Measurements of correlations between identical particles contain, besides amplitude in- 
terference, a plethora of other effects such as coherence, decaying resonances, variations 
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in impact parameter and momentum distribution, contamination by kaons and protons, 
final-state interactions etc. Understandably, a solid basis for subtracting all such ef- 
fects has been singularly hard to create. Given the number and degree of theoretical 
and experimental uncertainties entering even second-order correlations, the correspond- 
ing higher-order measurements have received only scant attention: if it is hard to extract 
source parameters in an honest and unambiguous way from second-order correlation data, 
it probably becomes even harder for third order. 

In the present paper, we elect to take a different approach: We measure higher-order 
correlations not so much with a view to extracting source parameters or "true" Bose- 
Einstein correlations, but in order to perform consistency checks. While few theorists 
have so far worked out the implications of their respective models for higher orders, this 
is in principle possible, and some examples of higher-order predictions exist || |3], |], ||, |6[ . 
If a given theory provides formulae for both second and higher orders, then these should 
apply to the corresponding data using the same parameter values throughout. 

With the aid of rapidly-improving measurement technology, we are attempting to put 
such predictions to the test. While the results reported here are quite preliminary in 
nature, they hopefully point the way to more general and sophisticated testing of theories 
rather than just measuring their respective parameters. Falsifying theories is arguably 
the best (some would say the only) way of making progress in a confused situation 0. 

Our tools for performing these consistency checks are cumulants and the correlation in- 
tegral |8], M. Cumulants, in subtracting out trivial lower-order contributions, have 
proven far more sensitive than the corresponding moments; their implementation in vari- 
ous forms of the correlation integral has, at the same time, improved statistical accuracy 
to a degree where such measurements have become meaningful. 



2 Quantum statistics theory 

The test we shall be reporting here is confined to one particular variable, the four- 
momentum difference = [{pi — pj) 2 — (Ei — Ej) 2 ] 1 ^ 2 . For this variable, the second 
and third-order cumulants are [TIJ 

C 2 (q) = p 2 (q) -pi®pi(g) , (1) 
C*3 = P3 ~ ^p2®pi + 2pi<8>pi<8>pi , (2) 

(3) 

where the third order quantities are functions of the three pair variables (912, ?23> fei)- 
These cumulants, including the crossed "<g>" quantities and event-mixing normalizations 
can be found from data samples in a precisely prescribed algorithm H]. 



The quantum statistics (QS) theory itself has a long and distinguished tradition [JT^, [TB| 



the version we concentrate on is based on analogies to quantum optics (for details, we 
refer the reader to Refs. [§, [| ||). Briefly, the main features of interest to us are: 
a) The pion field is split up into a "coherent" and a "chaotic" part: 

n(x) = n (x) + u ch { x ) . (3) 
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b) The ratio of chaotically created pions to the total number of pions is embodied in the 
"chaoticity parameter" , 

V = {rich) /{rich + n ) . (4) 

c) Much of the dynamics is contained within the normalized field correlator, 



(i4ft)iuS)> (i4ft-)iw%)> 



r _ \ L1 ch\ rv i) L1 ch^3)/ /p-x 

*J ~ r x .-> . .->... x .-> . .-> . -i 1/2 ' V / 



(where (A) = Tr(pA) is the ensemble average over states, weighted by the density matrix 
p) which is closely related to the Fourier transform of the chaotic field source functions, 
d) Working out two-point, three-point and higher-order averages, this theory of quantum 
statistics predicts unambiguously the normalized moments and cumulants of all orders. 
When relative phases are neglected, the first three "QS cumulants" of interest are [fj 

k 2 = -^- = 2p(l-p)d l2 +p 2 d 2 l2J (6) 

^3 = ~ — - — = 2p 2 (l - p)[d 12 d 23 + d 23 d 31 + d 31 d 12 ] + 2p 3 d 12 d 23 d 31 , (7) 
C 

k ± = ^ J ^ = JlP 3 ^ 1 ~P) d i2d 23 d 3i 

+ 2p A [di 2 d 23 d 3 idn + (^12^24^43^31 + ^14^42^23^3l] , (8) 

where the brackets under the sum indicate the number of permutations. These cumulants 
are functions of 1, 3 and 6 pair variables qij respectively. Note the combination of "ring"- 
and "snake" -like structures in the combinatorics. 

While in principle calculable from a given density matrix, the correlator is usually parametrized 
in a plausible and/or convenient way. Specifically, the parametrizations we shall be testing 
are, in terms of the 4-momentum difference correlators dij=d(qij), 

gaussian: d^ = exp(— r 2 q 2 -) , (9) 
exponential: d^ = exp(— rg^) , (10) 
power law: dy = q^ a . (11) 



3 UA1 data 

We have measured second- and third-order normalized cumulants using a sample of about 
160,000 minimum bias events taken with the UAl detector for pp collisions at 630 GeV/c. 
For details of the detector and other experimental information regarding particle pairs, the 



reader is referred to Ref. |14[ The following cuts were applied to this sample: —3 < 77 < 3, 
P± > 0.15 GeV, 45° < < 135° (by means of this "good azimuth" cut, our statistics were 
reduced considerably but acceptance corrections due to the "holes" in the UAl detector 
at small <ft were thereby avoided). Cumulants were calculated for positives and negatives 
separately and then averaged to yield "like-sign" values. No Coulomb corrections were 
applied. 
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Figure 1: Second order like-sign UA1 cumulant with fits using various parametrizations 
for d inside both the quantum statistical (QS) formula (§) and a simple parametrization 
AK 2 = pd. 

4 Fits to the second order cumulant 

In Figure 1, we show the second-order like-sign differential cumulant AK 2 = (/ p 2 j J Pi®P\) 
— 1, where numerator and denominator are integrals over bins spaced logarithmically 
between q = 1 GeV and 20 MeVQ. Fits to the data were performed using the three 
parametrizations (PD-flllD, either in the full QS form (P) or in a simple form AK 2 = pd\ 2 . 
All fits shown include, besides the free parameters p and r (or a), an additive constant as 
free parameter. These additive constants, necessary because UAl data is non-poissonian 
in nature, will be commented on further below. Best fit parameter values obtained were 
p = 0.66 ± 0.07, r = 1.16 ± 0.05 fm for the QS exponential and p = 0.05 ± 0.01, 
a = 0.64 ±0.05 for the QS power law. Goodness-of-fits were x 2 /NDF = 1.3, 4.2 and 11.5 
for QS power, exponential and gaussian respectively. 



To check its influence on fit values, the data point at smallest q, being of doubtful quality, 
was excluded; the resulting fit values do not differ much from the full fit. We note that 
the QS exponential misses the last three points (apart from the point at q = 20 MeV) and 
that the power laws (single or QS) appear to do the best job. The gaussian fits are too 
bad to warrant further attention and will be neglected from here on. Similar conclusions 
were reached by UAl earlier ||14j|. 



3 It should be remarked that previous work has shown the utility of using logarithmic rather than 
linear binning: much of what is interesting in correlations happens at small q, and this region is probed 
better by using logarithmic bins. 
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5 Consistency checks with third order cumulants 



As stressed already, we are interested not so much in obtaining numerical values for param- 
eters but rather in using these to check the theoretical formulae ©-(0) for consistency 
with the data. Three separate checks were performed: two based on approximations, the 
third involving a novel approach tentatively called "theory®experiment" which will be 
explained in Section 5.2. 



5.1 Approximate checks 

Third-order correlations and cumulants are functions of the three pair variables (qi 2 , ^23, 931), 
so that the question arises how best to view such three-dimensional correlations. The eas- 
iest projection involves setting the pair variables equal p, ^ [1^], qi2=q23=q3i=q, so that 
Eq. ([?D reduces to the simple formula 

k 3 (q) = 6p 2 (l - p)d 2 + 2p 3 d 3 . (12) 

Experimentally, however, the prescription of three mutually equal g's is so restrictive as to 
make measurement impossible. The usual way out ||T5| , |16|| has been to include all triplets 
whose mean of the three g^-'s is equal to a given q while still applying Eq. (|T2"| ) (the effect 
of this approximation has, to our knowledge, not been checked). 

The second approximation involves setting p = 1 without restricting the pair variables. 
Fortuitously, Eqs. (|7|), ([TDD then become k s = 2exp[— r(g 12 + ^23 + l3i)}, so that a simple 
change to the "GHP sum" variable S = (q , i2+Q , 23+?3i) does the trick. 

In Figure 2, we show the UAl third-order cumulant AK 3 as a function of the GHP sum 
variable S. The lower line represents the first approximation, i.e. formula (|12"D using 
the exponential parametrization ( |T0D and best-fit values from AK2 plus an arbitrary 
additive constant. (Similar approximations using the gaussian form (^) with the variable^ 

Q = \J q\ 2 + <?23 + ^li an d equal pair g's have been used before fll6fl .) The upper line, 
representing the second approximate check, was calculated by first fitting AK2 with p=l 
and an QS exponential for d to obtain r = 0.89 ±0.02 fm (not shown) and then importing 
this value into k 3 (p=l) = 2exp(—rS). 

We see that, in both cases, the theoretical curves lie well below the AK 3 data. Even 
an arbitrary shift by an additive constant does not improve the match because of the 
different shape of the curves as compared to the data points. 

5.2 The "theory® experiment" method 

The approximate consistency checks performed above are unsatisfactory for two reasons: 
first, because they rely on simplifications of the formulae which may be unwarranted, 
second, because they are suitable only for the exponential parametrization (or, using Q, 
for the gaussian equivalent). As shown above, however, the data for AK 2 , while not 

Note that the linear sum S is quite distinct from the Pythagorean sum variable Q. 
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Figure 2: Approximate test predictions, compared to the third-order UAl cumulant using 
GHP sum topology. 
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Figure 3: Third-order GHP max and GHP sum cumulants, together with theory®experi- 
ment predictions from QS theory and parameter values from AK2. Filled circles represent 
UAl data, triangles are predictions based on the QS power-law parametrization; squares 
are QS exponential predictions. 



6 



excluding an exponential form, would seem to prefer the power law — and the power law 
cannot be handled by these approximations. A seemingly better methodology emerges, 
surprisingly, from some considerations about normalization. 



Theory and theorists usually work with infinitesimally differential normalized quantities; 
for example, the second-order normalized cumulant is often written down as 

R(h, h) = P2i >H - 1 = C *&>fy. (13) 
Pi(h)pi(k 2 ) Pi(h)pi(k 2 ) 

which is (implicitly) fully differential in the momenta ki,k 2 . Similarly, the normalized 
theoretical cumulants ki = Ci/ pi® ■ ■ ■ ®p\ used above assume essentially perfect mea- 
surement accuracy and infinite statistics. 

Experimentally, one can never measure fully differential quantities; rather, the numerator 
and denominator are averaged over some bin of finite size Q (however small) before the 
ratio is taken; for example 

AK 2 (Sl) = /"^g , (14) 

which approaches the theoretical cumulant k 2 (q) only in the limit Q — > 0. 

This observation can be converted into an exact prescription for folding a given theoretical 
normalized quantity with experimentally measured one-particle distributions. For simplic- 
ity, we take second order quantities as an example. Since trivially C 2 (q) = k 2 (q) p\®p\(q), 
we can take k 2 from theory, p\®p\ from experiment and write exactly 



a v fo . In kf(q) Pi®P? P \q) dq _ / n Cf dq 

flA 2l"J - : — cxpt , : , = 7 „ expt/ \ j • 

In Pi® Pi (l) dq h Pi® Pi P (q) dq 



Correlation integral theory prescribes that || |10[ 

Pi®tF(q) = ((Y,8[Q-Q$\)), (16) 



where Qf- = [(pi 0, — p~j b ) 2 — (E? — Ej) 2 ] 1 ^ 2 is the four-momentum difference between two 
tracks i and j taken from different events a and b. Taking, for example, the QS cumulant 
(P) and the exponential parametrization (|l0|) , this leads to 

Cr expt (<?) = ((E% - Q^][Ml -p)eM~rQ%) +p 2 e W (~2rQ^ ^ , (17) 

which can be binned in q or otherwise integrated. In passing, we observe that Eq. flT5| ) 
reduces to the theoretical k 2 for infinitesimal Q or for constant pi®p\ as required. 

Clearly, this can be generalized to all possible moments and cumulants, independently of 
variable or integration topology. The procedure exemplified by Eq. ([17]) and its generaliza- 
tions amounts to a Monte Carlo integration of a theoretical correlation function sampled 
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according to the experimental uncorrelated one-particle distribution] for this reason, we 
like to call it by the diminutive "Monte Karli" or "MK" for short. MK can, of course, be 
implemented only for fixed numerical values of the theoretical parameters, in this case p 
and r. These must be determined either by more naive fitting methods (and then checked 
for consistency) or by a very cumbersome fitting procedure using the full event sample 
many times over. 

In Figure 3, the results of implementing the MK prescription are shown. Besides the 
GHP sum topology used in (b), we show in (a) a separate analysis using the "GHP 
max" topology ||10|| , which bins triplets according to the largest of the three variables, 
max(q , 12 , q 2 %, (fei)- Fit parameter values used for the respective power law and exponential 
MK points were taken from the naive QS fit to AK 2 of Figure 1. (The consistency of this 
procedure was checked by inserting these parameter values back into the MK formulae 
for AK 2 and finding agreement between UAl data and MK predictions.) Again, all MK 
points shown are determined only up to an additive constant, so that the curves may be 
shifted up and down. It is again clear, though, that the shape of third-order cumulant 
data measured differs appreciably from that predicted by the QS formulae and parameter 
values from AK 2 . This conclusion holds independently of the topology used and of the 
functional form taken for d. 



6 Discussion 

Concerning the fits to AK 2 , we have concluded that the gaussian parametrization dij = 
exp(—r 2 qfj) is quite unsuitable, while the exponential is better but not good. The best 
fit was obtained using either a simple or QS (double) power law. This confirms earlier 
results ||14||||17|| . The fits were reasonably stable even when excluding the point at smallest 
q, so that the effect is not due to this last point. 

Parameter values obtained from fits AK 2 were then applied to third-order cumulant data 
in three different checks. Both the two approximations as well as the exact theory®experiment 
(Monte Karli) prescription yielded predictions that did not match the data. The tests 
performed in this paper, namely checking three specific parametrizations (gaussian, ex- 
ponential and power-law) within one specific variable q for consistency between AK 2 and 
AK 3 appear to indicate that, under these specific conditions, the theory is contradicted 
by the data. 

It should be clear, though, that this conclusion can at this stage be preliminary and 
limited in scope only, for the following reasons: 

• The data shown is preliminary only and will have to await further checks such as 
acceptance corrections, full-azimuth studies, sensitivity to binning, etc. 

• The most important caveat relates to the structure of the overall multiplicity dis- 



tribution. The fact that UAl data is not poissonian in nature [[T8[ can be seen 
immediately at large q where AK 2 converges not to zero (as a poissonian cumulant 
would) but to ~ 0.4. The same holds for AK 3 . Theories, however, are almost uni- 
versally based on an overall poissonian: as can be easily verified from Eqs. 
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all cumulants tend to zero for large q. The policy followed here, namely reconciling 
poissonian theory with non-poissonian data by means of an additive constant in 
the cumulants, is a sensible but hitherto poorly-understood first step. The question 
of handling cumulants more adequately within a non-poissonian overall multiplic- 
ity distribution is presently being considered [jn|. We also hope that our results 
may goad theorists into more careful consideration of their work with respect to the 
implicit poissonian normalization used in most theories. See also Ref. . 



Closely related to these additive constants is the question of correct normalization. 
Traditional lore in second order divides the density correlation function (moment) 
AF 2 by an additional normalization factor /, taken as the moment at some large 
value of q. An alternative methods creates as many background pairs as necessary 
to achieve the limit of unity for AF 2 . While the third order moment can similarly 
be normalized to unity, the prescription fails for third order cumulants. A brief scan 
of the literature on third-order cumulants reveals that no adjustments were made 
for possible non-poissonian multiplicity structure. |3|, f|, |I| |9|, |K], [TT], [15], |2TJ, 



• Finally, one may mention possible changes to the present application of the theory 
such as inclusion of relative phases in the correlators, possible non-gaussian source 
currents, modelling the momentum dependence || of p, variable transformations 
p6|j in d and so on. 

Beyond these caveats, the following points are of relevance to the interpretation of our 
results: 

• No corrections for Coulomb repulsion |]J were included. We could argue that the 
same Coulomb effects that might shift AK 2 data upward would increase AK 3 data 
even more, since there are three pairs involved rather than one. Even more convinc- 
ing is the fact that AK% data rises more strongly than theoretical predictions even 
for large q (several hundred MeV) where Coulomb repulsion is not expected to be 
important. 

• Strictly speaking, the good power- law fit to AK 2 in itself is inconsistent: QS the- 
ory requires lim^o d(q) = 1, while the power-law parametrization diverges. 
Attempts to explain this in terms of variable transformations |2Bj or source size 



distributions p7|, g8| may therefore provide useful starting points in explaining the 
discrepancies in AK3. 

Track mismatching can lead to strong correlation effects because the reconstruction 
program may split a single track into a closely correlated pair. A great deal of effort 
in early experimental intermittency studies went into creating clean "split-track- 
finding" algorithms p9j and these are included in our analysis. We have checked 



through additional small-g cuts that mismatching does not appear to explain our 
strong rises in the cumulants. 

The UAl sample consists of ~ 15% kaons and protons which cannot be distinguished 
from pions. The effect that these would have on AK 3 is unclear. 
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• Resonances are known to increase AK2 at small q, the main effect in second order 
deriving from interference between "direct" pions and resonance-decay products.^ 
How and whether resonances would contribute to like-sign cumulants in third order 
(and for values of q of several hundred MeV shown in AK 3 ) is still quite mysterious. 

• At this point, one could wonder whether it is wise to even attempt to eliminate 
resonances from hadron-hadron collision data: apart from the theoretical and tech- 
nical difficulties, what dynamical information does the typical "size" r < 1 fm of 
a hypothetical "source" contain that is more important than a cascade structure 
containing resonances whose existence is beyond doubt? If the "source" is scarcely 
larger than a nucleon, then how can one speak of incoherent or even classical pro- 
duction of two pions? And if one eliminates long-lived resonances, then one would 
presumably still be left with the short-lived ones rather than the holy grail of an 
abstract quantum mechanical "source" . 

The results of the present paper may appear, at first sight, to contradict the conclusion 
||, based on an earlier UAl paper that QS theory was compatible with higher-order 



moments. The apparent discrepancy is explained by pointing out that 1) measurement 
techniques have improved considerably since then, 2) these techniques have permitted the 
present direct measurements of cumulants, which are considerably more sensitive than 
moments, and 3) even for these moment fits |15[], the radii were not quite constant but 
showed a systematic increase. 



Bose-Einstein correlation measurements with a view to extracting source parameters are 
by now well-established in hadronic and heavy ion phenomenology. Our intention here 
was to show that consistency checks between cumulants of different orders might be a 
second route to learning something about the system: if by this method a given theory 
can be tested already on a qualitative rather than quantitative basis, then opportunities 
for feedback and improvement of such theories may expand. 
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